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ABSTRACT 

Metal pollution in white dwarf atmospheres is likely to be a signature of remnant 
planetary systems. Most explanations for this pollution predict a sharp decrease in 
the number of polluted systems with white dwarf cooling age. Observations do not 
confirm this trend, and metal pollution in old (1-5 Gyr) white dwarfs is difficult to 
explain. We propose an alternative, time-independent mechanism to produce the white 
dwarf pollution. The orbit of a wide binary companion can be perturbed by Galactic 
tides, approaching close to the primary star for the first time after billions of years of 
evolution on the white dwarf branch. We show that such a close approach perturbs a 
planetary system orbiting the white dwarf, scattering planetesimals onto star-grazing 
orbits, in a manner that could pollute the white dwarf’s atmosphere. Our estimates 
find that this mechanism is likely to contribute to metal pollution, alongside other 
mechanisms, in up to a few percent of an observed sample of white dwarfs with wide 
binary companions, independent of white dwarf age. This age independence is the key 
difference between this wide binary mechanism and others mechanisms suggested in 
the literature to explain white dwarf pollution. Gurrent observational samples are not 
large enough to assess whether this mechanism makes a significant contribution to the 
population of polluted white dwarfs, for which better constraints on the wide binary 
population are required, such as those that will be obtained in the near future with 
Gaia. 

Key words: planet-star interactions, planets and satellites: dynamical evolution and 
stability, stars: evolution, stars: AGB and post-AGB, Oort Gloud, stars:kinematics 
and dynamics 


1 INTRODUCTION 

Elements heavier than helium sink rapidly in the atmo¬ 
spheres of white dwarfs, where sinking timescales are on the 
order of days to weeks (DA white dwarfs) and 10“* — 10® 
years (DB white dwarfs) (Koester & Wilken 2006) or see 
Fig. 1 of Wyatt et al. (2014). Observations of emission lines 
from metallic species in the atmospheres of white dwarfs 
(e.g. Koester et al. 1997; Zuckerman et al. 2003; Melis et al. 
2010), therefore, suggest the recent accretion of the ob¬ 
served material. Accretion from the interstellar medium 
was ruled out by Farihi et al. (2010); Aannestad et al. 
(1993); Jura (2006); Kilic & Redfield (2007); Barstow et al. 
(2014). These observations are thought to indicate the 
accretion of planetary material (e.g. Alcock et al. 1986; 
Debes & Sigurdsson 2002; Jura 2003, 2008; Kilic et al. 
2006), in at least 25% (Zuckerman et al. 2003, 2010), but 
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maybe up to 50% (Koester et al. 2014) of white dwarfs. 
Dusty and/or gaseous discs, very close (r < Rq) to a subset 
of these polluted white dwarfs may be indicative of the ac¬ 
creting material (e.g Kilic et al. 2006; Gansicke et al. 2006; 
von Hippel et al. 2007; Jura et al. 2007; Gansicke et al. 
2007). The composition of the accreted material, when anal¬ 
ysed in detail supports the planetary material hypothesis 
(e.g. Klein et al. 2010, 2011; Gansicke et al. 2012). 

Planetary systems are common on the main-sequence, 
and throughout the Milky Way they are the rule, not 
the exception (Cassan et al. 2012). Observations with 
HARPS suggest that at least 50% of solar-type stars har¬ 
bour at least one planet with a period less than 100 
days (Mayor et al. 2011), whilst Herschel found that at 
least 20% of FGK stars have a detectable debris disc 
(Eiroa et al. 2013). Although close-in (< 1 — 5AU) plan¬ 
ets may be swallowed by the expanding giant’s enve¬ 
lope (Villaver & Livio 2007, 2009; Mustill & Villaver 2012; 
Adams & Bloch 2013; Villaver et al. 2014), there is no ev- 
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idence that all outer planetary systems are destroyed (e.g. 
Jura 2004; Bonsor & Wyatt 2010). Dynamical instabilities 
in the outer planetary system, induced following stellar 
mass loss (e.g Debes & Sigurdsson 2002; Veras et al. 2011; 
Veras & Wyatt 2012; Veras & Tout 2012; Veras et al. 2013c; 
Voyatzis et al. 2013; Mustill et al. 2014; Veras & Gansicke 
2015), may scatter asteroids or comets onto star-grazing or¬ 
bits (Bonsor et al. 2011; Bonsor & Wyatt 2012; Debes et al. 
2012; Frewen & Hansen 2014), where they are tidally dis¬ 
rupted and accreted onto the star (e.g. Graham et al. 1990; 
Jura 2003, 2008; Bear & Soker 2013; Veras et al. 2014c, 
2015c). 

Given a sufficiently massive planetesimal belt, it has 
been shown that such a theory can explain the ob¬ 
served metal pollution (Bonsor et al. 2011; Debes et al. 
2012; Frewen & Hansen 2014). However, a steep decrease 
in the mass of asteroids/comets scattered onto star-grazing 
orbits with time after the formation of the white dwarf, is 
predicted (Bonsor et al. 2011; Debes et al. 2012). This find¬ 
ing would correspond to a steep decrease in the level of pol¬ 
lution (or fraction of systems with pollution), as a function 
of time. Although late time instabilities may be induced in 
multi-planet systems (Veras et al. 2013c; Mustill et al. 2014; 
Veras & Gansicke 2015), the frequency of these will always 
be less than those around young white dwarfs. 

Observationally there is no evidence that there are 
fewer, old polluted white dwarfs, nor that the level of pollu¬ 
tion decreases with white dwarf cooling age (Wyatt et al. 
2014; Koester et al. 2014). Even the archetypal polluted 
white dwarf, van Maanen’s star (van Maanen 1920), has 
an effective temperature of 6,220K, or cooling age of ~ 
3Gyr (Sion et al. 2009), and there are many further ob¬ 
servations of polluted, old, white dwarfs (Teff < 8, OOOK, 
equivalent to cooling ages of up to 5 Gyrs Farihi et al. 
(2011); Koester et al. (2011)). Instabilities induced follow¬ 
ing stellar mass loss struggle to explain pollution for white 
dwarfs with cooling ages of Gyrs (Debes & Sigurdsson 2002; 
Bonsor et al. 2011; Debes et al. 2012). Alternative explana¬ 
tions for the pollution in these systems have been suggested, 
including stellar encounters (Farihi et al. 2011), or a rela¬ 
tionship with magnetic helds (Hollands et al. 2015). Others, 
such as exo-Oort cloud comet impacts (Veras et al. 2014d; 
Stone et al. 2015) and volatile sublimation of minor bodies 
in planet-less systems (Veras et al. 2015a), have largely been 
ruled out. 

Here, we suggest an alternative explanation. It has been 
shown for main-sequence planetary systems, that whilst in 
general wide (a > 1,000AU) binary companions do not in¬ 
fluence the dynamics of the planetary system, the orbit of 
the binary varies due to Galactic tides. During periods of 
‘close’ pericentre passage, the binary may excite the eccen¬ 
tricities of planets, even ejecting them (Kaib et al. 2013). 
The increased eccentricity of exoplanets in systems with 
wide binary companions has been confirmed by observations 
(Kaib et al. 2013). Here, we propose that the same scenario 
could be applied to white dwarf planetary systems. If the 
white dwarf is orbited by a wide binary companion, the plan¬ 
etary system may remain unperturbed for billions of years, 
before Galactic tides alter the orbit of the companion such 
that it induces dynamical instabilities in the planetary sys¬ 
tem. These instabilities can lead to material being scattered 
onto star-grazing orbits and accreted onto the white dwarf. 


such that we observe pollution in the white dwarf’s atmo¬ 
sphere. This process has the advantage of being independent 
of the age of the white dwarf, although it depends critically 
on the population of wide binaries orbiting white dwarfs. 

The population of wide (> 1,000 AU) binaries 

is difficult to characterise, in general. Common proper 
motions indicate that two stars may be related {e.g. 
Lepine & Bongiorno 2007; Makarov et al. 2008), but follow¬ 
up parallax and radial velocity observations are required 
to determine whether they are bound. Wide binaries are, 
therefore, likely to be significantly more common than 
determined observationally. Theoretical models predict a 
wide binary fraction of 1-30%, with 10® < a < O.lpc 
(Kouwenhoven et al. 2010), whilst observations of solar-type 
stars (Lepine & Bongiorno 2007; Raghavan et al. 2010) find 
a fraction of 9.5% or 11.5%, and observations of field white 
dwarfs yield a wide binary fraction of 22% (Farihi et al. 
2005). There are, however, very few white dwarfs with known 
wide binary companions that have been searched for pollu¬ 
tion, and many surveys for pollution have focussed on ap¬ 
parently single white dwarfs. Zuckerman (2014) present a 
sample of 17 white dwarfs with companions separated by 
more than 1,000AU, searched for Ca H with Keck/VLT 
and find that 5 are polluted. We note here that white 
dwarfs in close binaries may be polluted by interactions 
with the stellar wind of the companion, as discussed in {e.g. 
Zuckerman et al. 2003; Zuckerman 2014), and for this reason 
we focus on ‘wide’ binaries. Fig. 8 of Veras & Tout (2012) 
shows that the boundary between ‘close’ and ‘wide’ occurs 
at tens of AU. 

In this paper, we illustrate the manner in which a wide 
binary companion could lead to pollution in a white dwarf. 
This mechansim is likely to be just one of many that lead 
to pollution. The critical difference of this mechanism is the 
absence of any dependence on white dwarf cooling age. §2 
shows how a wide binary’s orbit can be altered by Galactic 
tides, such that a close approach occurs between the sec¬ 
ondary and a planetary system orbiting the primary (white 
dwarf). §3 presents simulations that show how this close ap¬ 
proach can lead to planetesimals scattered onto star-grazing 
orbits, a pre-requisite for pollution. In §4 we estimate the 
fraction of white dwarfs with wide binary companions where 
pollution might occur, which is compared to observations in 
§5. In §6 and §7 we discuss our results and present our con¬ 
clusions. 


2 THE CHANGE IN ORBIT OF A WIDE 
BINARY DUE TO THE GALACTIC TIDE 

We firstly illustrate the manner in which the Galactic tide 
can change the orbit of a wide binary. Above a critical sep¬ 
aration, the Galactic field strength is different for the two 
components of the binary, such that their orbit evolves. An¬ 
alytically the effect of the Galactic tide on the binary’s orbit 
can be understood through the equations of motion for the 
binary’s orbital parameters. Although we are unaware of a 
complete analytic solution to these equations (1-4), numer¬ 
ical solutions can be used to explore the phase space and 
help us to understand the evolution of the binary’s orbit. 
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2.1 Equations of motion 

As the majority of observed polluted white dwarfs are nearby 
{e.g. within O.lkpc Holberg et al. 2008), we restrict our¬ 
selves to the Solar neighbourhood, which we assume re¬ 
sides within 8kpc of the Galactic centre. This restriction 
allows us to neglect planar Galactic tides (see Fig. 2 of 
Veras & Evans 2013a), which represents a widely-used sim¬ 
plification (Heisler & Tremaiire 1986; Matese & Whitman 
1989, 1992; Matese et al. 1995; Breiter et al. 1996; Brasser 
2001; Breiter & Ratajczak 2005). Further, we coirsider only 
bound binary systems. Boundedness effectively equates to 
adiabaticity in the Galactic tidal regime because no escape 
can occur in the adibatic regime (Veras & Evans 2013a) and 
the non-adibatic limit is close to the system’s Hill ellipsoid 
(Veras et al. 2014a). For very wide binaries (at >> 10®AU) 
adiabicity is no longer a good approximation (shown by Fig. 
3 of Veras & Evans (2013b)) and the full equations of mo¬ 
tion describing this regime are shown in Column 5 of Table 
1 of Veras et al. (2014a). 

Here, we consider the adiabatic regime. Consequently 
we can, by considering vertical tides alone and using the 
perturbed two-body problem, derive the equations of mo¬ 
tion for the orbital elements of the binary (as shown in 
Veras & Evans (2013a)); semi-major axis, at, mean motion, 
Ub, eccentricity, et, inclination of the binary orbital plane to 
the Galactic plane, ib, argument of pericentre, uJb : 


dub 

dt 





( 1 ) 


5eby/l -e^ . . 2 . .Y- 

-^ cos sin (Uf, sm ibTzz (2) 




5e^ sin 2ujb sin 2ib 
5 sin^ UJb (sin^ ib — el) — (l — el) 

Vi^l 


T, 


(3) 


(4) 


The variable Tjj is determined by the Galactic model 
used, and contains information about the matter density and 
gravitational potential. Here, we use the disc component of 
the three-component model of Veras et al. (2013a). Because 
we fix the distance from the Galactic centre at 8 kpc, we 
also fix = -5.352 x 10“®“s“^ 


2.2 An example binary 

Although we are unaware of a complete analytic solution to 
equations (1-4), we can obtain numerical solutions that en¬ 
able us to follow the evolution of example orbits. There are 
several key properties of the evolution of binary orbits under 
the influence of the Galactic tide. Of particular relevance to 
this work is the evolution of the binary’s semi-major axis and 
eccentricity. The semi-major axis remains constant (Eq. 1), 
under the adiabatic approximation of these equations. The 
eccentricity evolves periodically, with a time period and am¬ 
plitude that depend on the initial conditions of the binary’s 
orbit, as well as the strength of the Galactic tide (T^j). If a 
binary’s orbit is to perturb the primary’s planetary system, 
we are interested in the maximum eccentricity that its orbit 
reaches and the times between eccentricity peaks. Thus, we 
focus on these two parameters in the following sections. 


Firstly, in order to illustrate our proposed scenario, let 
us consider, for example, a IMq binary companion orbiting 
a IMq primary, in the solar neighbourhood {R = 8kpc), 
with a semi-major axis of 3000AU, inclined to the Galac¬ 
tic plane by 80°. We gave the binary an initial eccentricity 
of eo = 0.2, and a range of initial arguments of pericen¬ 
tre (cuo = 5,40,90,160°) and followed the evolution of the 
binary’s eccentricity, shown in the top panel Fig. 1. The ec¬ 
centricity of this binary’s orbit evolves on time periods that 
are significantly longer than the main-sequence lifetimes of 
the stars, or even the age of the Universe. Therefore, if the 
binary starts at low eccentricity, as we assume, it remains 
at low eccentricity throughout the primary’s main-sequence 
lifetime. 

We then consider the fate of this binary system, as the 
primary becomes a white dwarf. So far, the primary’s plan¬ 
etary system is relatively unperturbed by the binary com¬ 
panion. By using a typical mass loss prescription^, the pri¬ 
mary becomes a white dwarf of mass 0.52Mq (see caption of 
Fig.4 of Veras et al. (2013c)). At this semi-major axis mass 
loss is likely to be non-adiabatic (see for example Fig.3 of 
(Veras et al. 2014a)), thus the exact expansion of the bi¬ 
nary’s orbit is difficult to predict (Veras et al. 2011). We, 
therefore, make a reasonable assumption that the binary’s 
orbit expands by a factor of 3 to at = 9,000AU. Both the re¬ 
duction in mass of the binary and the increase in semi-major 
axis of the binary, increase the ability of Galactic tides to 
change the binary’s orbit. As is shown in the bottom panel 
of Fig. 1, the binary’s eccentricity now evolves on ‘shorter’ 
timescales, several rather than tens of Gyrs. Thus, with our 
assumption that the binary’s orbit starts the white dwarf 
phase with an eccentricity of eo = 0.2, depending on cjq, it 
evolves to its maximum eccentricity (minimum pericentre) 
within 3-7Gyrs. Approaching the primary at a distance of 
closest approach of around 150AU, the binary is likely to 
perturb the primary’s planetary system. It is these pertur¬ 
bations that we invoke for their ability to produce polluted 
white dwarfs, particularly at late times. 


2.2.1 Including mass loss 

Fig. 1 integrates Eqs. 1-4 for the main-sequence and the 
white dwarf phases separately. It is possible to integrate 
these equations self-consistently, if a formulation for stellar 
mass loss is included. However, it is computationally inten¬ 
sive and therefore, we only perform this calculation for one 
example system. Using the stellar evolution prescription of 
Hurley et al. (2000), coupled with N-body integrations, as 
in Veras et al. (2011), Fig. 2 shows the evolution of the peri¬ 
centre of a binary orbit that starts with a 2Mq and a IMq 
planet on an orbit with Ob = 4, OOOAU, but where the 2Mq 
primary evolves to become a white dwarf of mass O.64M0. 
The same behaviour as in Fig. 1 is seen. The binary’s orbit 


^ The SSE code (Hurley et al. 2000) applies Reimer’s mass loss 
(Kudritzki & Reimers 1978) at early evolutionary stages, but 
during the AGB it applies the semi-empirical mass-loss rate of 
Vassiliadis &; Wood (1993), reaching a maximum during the su¬ 
perwind of 1.36 X lO®(L*/L0)MQyr“^, where L* is the stellar 
luminosity. 
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Figure 1. The change in orbital parameters due to Galactic 
tides and stellar evolution, for several example binaries. The top 
panel represents the primary on the main-sequence and the bot¬ 
tom panel, once the primary has evolved to become a white dwarf. 
Tides can produce extreme eccentricity oscillations during the 
white dwarf phase which do not occur along the main sequence, 
thereby triggering planetary scattering which can ultimately lead 
to white dwarf pollution. Assumed in these plots is that the bi¬ 
nary orbit expands non-adibatically by a factor of 3 due to mass 
loss from the secondary star M 2 . 

is not perturbed significantly during the main-sequence evo¬ 
lution, but undergoes long period oscillations once the star 
becomes a white dwarf. 

The mass loss is modelled to be isotropic, which rep¬ 
resents an excellent approximation to reality (Veras et al. 
2013b). In this case, the pericentre can never decrease dur¬ 
ing giant branch mass loss (Eq. 21 of Veras et al. (2011)). 


3 THE POLLUTION OF THE WHITE DWARF 

We have shown that Galactic tides can perturb orbits such 
that a secondary companion approaches close to any plane¬ 
tary system orbiting the primary (white dwarf). In this sec¬ 
tion we present simulations that illustrate how the close ap¬ 
proach of the secondary could perturb the planetary system 
orbiting the white dwarf in such a manner that planetesi- 
mals are scattered onto star-grazing orbits. For the purpose 
of this work we consider such scattering sufficient to produce 
white dwarf pollution, see §6 for a discussion of this assump¬ 
tion. These simulations are intended as a proof of concept, 
rather than a detailed analysis. 



Figure 2. A demonstration of how the first close approach in a 
binary system can occur during the white dwarf phase. Plotted 
is the self-consistent pericentre distance evolution of one partic¬ 
ular binary system through all phases of stellar evolution. The 
a initial binary parameters are: Mi = 2Mq, M 2 = IMq, ei, = 0.1, 

af, = 4000 au, if, = 110°, uif, = 0°. The red, blue and purple 
curves represent the main sequence, giant branch, and white dwarf 
phases of evolution, respectively. Adiabatic Galactic tidal equa¬ 
tions are used for the main sequence and white dwarf phases; the 
full non-adiabatic equations for both tides and mass loss are used 
along giant branch phases, with /o = 150°. 


3.1 Details of our simulations 

Clearly, the parameter space available to investigate, both 
in terms of the diversity of planetary system architectures 
and binary orbits, is enormous. We can only investigate ex¬ 
ample systems. We, therefore, consider a binary in which 
a IM 0 primary is orbited by a companion of mass Q.8Mq. 
The IM 0 primary evolves to become a O.52M0 white dwarf 
(Veras et al. 2013c). O.SMq represents a typical value for 
a main-sequence stellar mass (Parravano et al. 2011), and 
such stars will remain on the main-sequence for the current 
age of the Universe (~ 14 Gyrs). 

The primary is orbited by a very simplistic planetary 
system with one planet and one planetesimal belt or de¬ 
bris disc, similar to Neptune and the Kuiper belt in our So¬ 
lar System. The presence of both the planet and planetes- 
imals means that perturbations to the orbit of the planet 
can lead to planetesimals being scattered onto star-grazing 
orbits (Bonsor et al. 2011; Debes et al. 2012). In addition, 
the binary companion can perturb planetesimals directly 
onto star-grazing orbits (e.g. Marzari et al. 2005). We fix 
the plane of the planetary and binary orbits to be the same. 
We fix the planet’s semi-major axis at 30AU, which mim¬ 
ics Neptune’s semi-major axis, and the planetesimal belt 
runs from 30-50AU. We assume that at such a small orbital 
distance stellar mass loss is adiabatic, such that the white 
dwarf is orbited by a planet at 60AU and belt running from 
60-100AU. We neglect the potentially destructive effects of 
YORP spin up on the planetesimals (Veras et al. 2014b) and 
the potentially significant movement of these asteroids due 
to the Yarkovsky effect from giant branch stellar evolution 
(Veras et al. 2015b). These effects would have a greater in¬ 
fluence on an exo-asteroid belt located within lOAU than an 
exo-Kuiper belt located at several tens of AU. 

Our simulations are run using the Mercury N-body in- 
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tegrator, RADAU (Chambers 1999). We track particles that 
are scattered close to the star. It is these particles that we 
assume are tidally disrupted and pollute the white dwarf. 
The RADAU integrator is used in order to track the evolu¬ 
tion of these particles to small pericentre, in a manner that 
the sympletic integrator cannot. However, in order to limit 
the runtime for our simulations, we only track particles as 
close in as 0.1 AU. We consider the particles that are scat¬ 
tered interior to 0.1 AU to be representative of the popula¬ 
tion scattered onto star-grazing orbits where they would be 
disrupted, more likely to occur when r Rq (Veras et al. 
2014c). Given that our purpose is to illustrate the feasibility 
of this mechanism, we consider this assumption to be suf¬ 
ficient. Our choice merely means that we over-estimate the 
population scattered onto star-grazing orbits, but do so in a 
consistent manner in all our simulations, thereby removing 
comparison inconsistencies between simulations. 

We use our simulations as a tool to investigate the bi¬ 
nary orbits that perturb the planetary system and scatter 
particles onto star-grazing orbits. Therefore, we focus on 
short integrations (10-100 binary orbits) for which the bi¬ 
nary orbit is fixed and does not evolve due to the Galactic 
tide. These integrations can represent the evolution of the 
planetary system directly following the binary’s evolution 
onto the orbit considered. 

The exact orbits of the planet and planeteslmals will 
clearly evolve with time and the evolution of the binary’s 
orbit. It is not possible to consider all possible evolutionary 
paths. Instead we investigate the increase in the scattering 
of planetesimals onto star-grazing orbits in comparison with 
a simulation in which the binary is not present. Although 
this comparison may not always mimic the potential effects 
of an excited planetesimal belt, or scattered disc, it illus¬ 
trates clearly the perturbative effect of the binary. We start 
our planetesimals with orbits that are initially unperturbed 
by the planet and have semi-major axes, eccentricities, in¬ 
clinations, longitudes of pericentre, longitudes of ascending 
nodes and mean anomalies, which are randomly drawn from 
uniform distributions with ranges of 60 < a < lOOAU, 
0 < e < 0.02, 0 < i < 1° and 0° < A, D, M < 360°. 

3.2 An example system 

The aim of this section is to illustrate the manner in which 
the change in orbit of a binary due to Galactic tides can lead 
to the pollution of a white dwarf, even a white dwarf with 
a cooling age of Gyrs or longer. We consider the example 
binary system shown in Fig 1. The binary companion orbits 
at fli, = 3, OOOAU and whilst it is on the main-sequence has 
its eccentricity perturbed very little by the Galactic tide, 
up to an absolute maximum of 66 = 0.6, from 66 = 0.2 ini¬ 
tially, during a long main-sequence lifetime of lOGyr. Even 
if the binary remains on an orbit with 66 = 0.6 for 900Myr, 
a typical main-sequence lifetime, an outer belt survives or¬ 
biting the primary star. In fact, our test simulation finds 
that more than 60% of the original belt remain on stable or¬ 
bits. Given that Fig. 1 shows that the binary’s orbit evolves 
and will only be on such an eccentric orbit for timescales 
significantly shorter than the primary’s full main-sequence 
lifetime. Thus, in this example, we show that a planetary 
system can survive the star’s main-sequence lifetime with¬ 
out being perturbed unduly by its binary companion. 


1.01 . 
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Figure 3. The semi-major axis, eccentricity distribution of plan¬ 
etesimals orbiting the white dwarf, prior and following the peri¬ 
centre passage of a wide binary companion with 66 = 0.99. This 
is the same system as shown in the bottom panel of Fig. 1. This 
figure demonstrates that despite being dynamically excited, an 
exo-Kuiper belt largely survives stellar evolution of a primary 
star with a wide-orbit companion. 


Once the primary star in this example binary evolves to 
become a white dwarf, the orbits of the binary and the pri¬ 
mary’s planetary system expand, for the planetary system 
adiabatically, but for the binary non-adiabatically. As above, 
we estimate that the binary’s orbit moves to as = 9, OOOAU 
and consider the perturbations that it now makes on the pri¬ 
mary’s planetary system, that has expanded adiabatically by 
a factor of 2 (belt is now at about 60-100AU). The binary’s 
orbit is now perturbed by the Galactic tides, evolving up to 
an eccentricity of 0.978 (q ~ 200AU) after several Gyrs of 
evolution on the white dwarf branch. As the binary’s eccen¬ 
tricity increases, so does its influence on the planetary sys¬ 
tem, particularly during pericentre passages. Fig. 3 shows 
the effect on our example planetary system of a binary with 
eccentricity, 66 = 0.99. Fig. 3 shows that prior to the binary’s 
pericentre passage, the planetary system is relatively undis¬ 
turbed, but that following its pericentre passage the planet 
and many planetesimals have been scattered. As discussed 
above, the scattered planetesimals have the potential to pol¬ 
lute a white dwarf. Fig. 4 shows the number of planetesimals 
scattered onto star-grazing orbits following a pericentre pas- 
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Figure 4. Our example white dwarf planetary system, orbited 
by a binary companion with = 9, OOOAU (shown in Fig. 1 
and Fig. 3). The number of particles (N) that hit the star, as a 
function of time, in units of the binary’s orbit period, with t = 0, 
occurring at pericentre. The high binary eccentricity induced by 
the Galactic tide increase the number of planetesimal collisions 
with the white dwarf. 

sage by a binary on orbits with various eccentricities. For 
binary eccentricities above 0.8, significantly more planetesi- 
mals are scattered onto star-grazing orbits than without the 
presence of the binary. If we refer back to Fig. 1, depending 
on the exact parameters of the binary’s orbit, the binary’s 
eccentricity could be increased from 0.2 to above 0.8 be¬ 
tween 2 and 6 Gyrs after the start of the white dwarf phase. 
In this manner we illustrate a mechanism that could lead to 
the pollution of old white dwarfs. 

3.3 A brief investigation of the parameter space 

The above section illustrated one example system in which 
white dwarf pollution could occur in white dwarfs with cool¬ 
ing ages of Gyrs. Clearly, the parameter space of binary or¬ 
bits is huge and this mechanism will not work for all of them. 
A detailed investigation of the entire parameter space is be¬ 
yond the scope of this work, and not necessary to illustrate 
the potential of this mechanism. Instead, we perform a small 
additional suite of N-body simulations in order to show that 
there are many binary orbits for which this mechanism could 
work. We focus here on the white dwarf phase only. 

We consider binary companions with semi-major axes 
of at = 2, 500AU, at = 5, OOOAU and at = 9, OOOAU. The 
eccentricity of the binary companion is varied, but fixed 
throughout each individual simulation. Fig. 5 shows the 
number of planetesimals scattered onto star-grazing orbits 
as a function of time. The pericentre passages of the binary 
are evident from the increase in the scattering rate, creating 
sawtooth-like profiles. As noted in §3 these numbers must 
be compared to the simulation with the binary removed, 
in order to remove some dependencies on the initial con¬ 
ditions used. For sufficiently high eccentricity companions, 
it is clear that more material hits the sun than without a 
binary companion. A conservative approximation to the ec¬ 
centricity at which this occurs can be obtained from the 
criterion described in Holman & Wiegert (1999) for the sta- 
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Figure 5. The number of particles (N) scattered onto star¬ 
grazing orbits as a function of time, for at = 2,500AU (top), 
at = 5, OOOAU (middle) and at = 9, OOOAU (bottom). The planet 
is ejected tor et = 0.99, at = 9, OOOAU after O.SMyr, hence the 
scattering stops. This figure demonstrates that as the binary’s 
orbit evolves to high eccentricity, the rate at which particles are 
scattered inwards increases significantly. 

bility of planetary systems orbiting binary systems. Fixing 
acrit at 60AU in their Eq.l, although technically only valid 
for et < 0.8, yields approximately the same critical eccen¬ 
tricities for the binary (et > 0.8; at = 2,500AU; et > 0.9, 
at = 5, OOOAU; et > 0.95,06 = 9, OOOAU) as found from our 
simulations (see Fig. 6). Alternative criteria exist for the 
stability of three body systems (Eggleton & Kiseleva 1995; 
Petrovich 2015), and although valid for higher eccentricity, 
have not been derived for two planets orbiting one star. 
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Figure 6. The analytic criterion (solid line) for the stability of cir- 
cumprimary orbits from Holman & Wiegert (1999) (Eq. 1). The 
arrows show a comparison with the results of our simulations for 
the white dwarf phase, on short timescales. The shaded region in¬ 
dicates the limit of validity of this approximation, e > 0.8, which 
unfortunately corresponds to the region of interest. 

rather than a binary companion and a planet orbiting the 
central star. 

The critical point to take from these simulations is that 
planetesimals are scattered onto star-grazing orbits whilst a 
binary companion has high eccentricity, but remain on more 
stable orbits whilst the binary’s eccentricity is lower. Thus, 
planetesimals may survive the star’s main-sequence evolu¬ 
tion, as well as some proportion of its white dwarf evolution 
on stable orbits, before being scattered onto star-grazing or¬ 
bits as the binary’s orbit evolves to high eccentricity. 

3.4 Caveats and Assumptions 

The purpose of the above simulations is to illustrate that a 
wide binary companion could lead to white dwarf pollution 
at late times by scattering planetesimals onto star-grazing 
orbits. This requires that sufficient material survives to the 
white dwarf phase, as well sufficient material is scattered 
onto star-grazing orbits to produce detectable pollution. 

We assume here that planetesimals scattered onto star¬ 
grazing orbits lead to white dwarf pollution. This is a 
fairly robust assumption, as discussed in previous work, 
although gaps still exist in our knowledge of the exact 
processes involved (e.g. Jura 2003; Bonsor et al. 2011; 
Debes et al. 2012; Frewen & Hansen 2014). Planetesimals 
scattered close to the star are tidally disrupted, and ac¬ 
crete onto the star (Veras et al. 2014c, 2015c), potentially 
via dusty or gaseous accretion discs (Hartmann et al. 2011; 
Rafikov 2011a,b; Bochkarev & Rafikov 2011; Metzger et al. 
2012; Rafikov & Garmilla 2012). Our simulations are good 
at showing that material is scattered onto star-grazing or¬ 
bits, but only indicate increased rates of scattering, rather 
than exact levels, and cannot, therefore, be used to predict 
accretion rates. However, the amount of material required 
to produce even the most polluted systems (e.g. a Ceres 
mass for Dufour et al. (2010) or see Fig. 9 of 

Girven et al. (2012)) is low, and our simulations indicate 
high levels of scattering (Fig. 5). These masses are low com¬ 


pared to typical planetesimal belt masses thought to ex¬ 
ist on the main-sequence, even those below detection limits 
(Wyatt et al. 2007; Wyatt 2008), such as the Kuiper belt 
with its mass of around O.IM 0 (Gladman et al. 2001). Al¬ 
though these masses may be reduced by collisions, or dynam¬ 
ics as material is scattered during the star’s main-sequence 
lifetime, we do not anticipate that in most systems these 
processes remove sufficient material to prevent detectable 
pollution being produced (Bonsor & Wyatt 2010), although 
of course in some planetary systems dynamical clearing may 
be highly efficient {e.g. Raymond et al. 2009, 2010). 

We choose to only simulate binaries with fixed eccentric¬ 
ity, in order to simplify our results and save computational 
time. From this choice we can clearly show an increased rate 
of scattering relative to simulations where the binary is not 
present. However, we miss any changes due to the exact evo¬ 
lution of the system as the binary’s orbit increases smoothly 
in eccentricity, or due to particularities of the exact configu¬ 
ration of the system {e.g. planetesimals with excited eccen¬ 
tricities or inclinations or trapped in resonance etc ). These 
details may change the scattering rates. However, given our 
vague knowledge of the structure of any planetary system, 
we consider that the present simulations are sufficient to 
show that material is scattered and a white dwarf could be 
polluted in this manner. 


4 ESTIMATING THE FREQUENCY OF 

POLLUTION CAUSED BY A WIDE BINARY 

COMPANION 

Having shown that for particular orbital parameters of the 
binary and planetary system, pollution of the white dwarf 
can occur, it would be useful to estimate the contribution of 
this mechanism to the population of polluted white dwarfs. 
This value firstly depends on the binary having an appropri¬ 
ate orbit that remains at low eccentricity for the primary’s 
main-sequence evolution, but increases to sufficiently high 
eccentricity to perturb the primary’s planetary system dur¬ 
ing the primary’s white dwarf evolution, and secondly, on 
the primary having a planetary system with sufficient ma¬ 
terial that can be scattered onto star-grazing orbits in order 
to produce the pollution. 

4.1 The binary orbit 

Firstly, we note that the binary orbits that are most likely 
to be interesting in terms of the proposed scenario are those 
at intermediate wide orbits {i.e. 3,000 ^ as < 10,000AU). 
Here, the influence of Galactic tides is strong, but not suf¬ 
ficiently strong that the main-sequence planetary system 
would have been perturbed (as might occur for larger semi¬ 
major axes) and the binary does not orbit too close to the 
primary, so that the primary’s planetary system is only per¬ 
turbed when the binary’s orbit increases significantly in ec¬ 
centricity. 

Here, we estimate the fraction, /wb, of wide binaries 
where the eccentricity of the binary’s orbit increases suf¬ 
ficiently (above a critical value), during the white dwarf 
phase, such that the primary’s planetary system is per¬ 
turbed. The distribution of orbital parameters of wide bina¬ 
ries in the galaxy are not well known. We, therefore, con- 
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Figure 7. Estimating the fraction of systems where the wide bi¬ 
nary mechanism has the potential to pollute the primary’s plan¬ 
etary system. Plotted is the cumulative fraction of systems that 
have a pericentre passage closer than a critical value during the 
white dwarf phase. The maximum eccentricity (minimum peri¬ 
centre) can be estimated from Fig. 5; see discussion in §4.1 for 
details. 

sider a uniform distribution of initial binary eccentricity, 
eo , longitude of pericentre, uq, sine of inclination to the 
Galactic plane, sin to, and use a Monte-Carlo model to se¬ 
lect 300 binary orbits for each semi-major axis^. We fix the 
binary semi-major axis at 06 = 9, OOOAU, aj, = 2, 500AU or 
06 = 5, OOOAU and the masses of the pair fixed at Mwd = 
0.52Mq (evolved from mi = IMq ) and M 2 ~ O.SMq), 
in line with Fig. 1 and Fig. 5. Fig. 7 shows the fraction 
of systems where a pericentre passage closer than a critical 
value (or equivalently, the eccentricity evolves to above a 
critical value), at any point during the white dwarf phase. 
A maximum time period of 14Gyr, the age of the Universe, 
is considered. From Fig. 5, an eccentricity of greater than 
approximately 0.8 {ah = 2,500AU), 0.9(a6 = 5,OOOAU) or 
0.95 (ah = 9, OOOAU) is required to perturb a planetary sys¬ 
tem at 60-100AU. In other words, the binary must enter the 
region interior to ~ 500AU to perturb the planetary system. 
Fig. 7 shows that 35%, 17% and 10% of systems evolve to 
this pericentre. Although clearly this value will vary signif¬ 
icantly with the binary’s semi-major axis, we consider 20% 
to be a reasonable estimate for the fraction of wide bina¬ 
ries with orbital parameters in the correct range to produce 
polluted white dwarfs. 

Another constraint on the orbit of the binary results 
from the requirement that a planetary system survives, rel¬ 
atively unperturbed, orbiting the primary until late times. 
This requirement means that the binary’s eccentricity must 
remain below a critical value during the primary’s main- 
sequence evolution. This value can be roughly estimated 
by considering the criterion of Holman & Wiegert (1999) 
(Eq 1.) for stable circumprimary orbits. For example for 
our three example binaries with ah = 2, 500AU, 5,OOOAU 
and 9,OOOAU, the binary eccentricity must remain below 
66 < 0.4, 0.67 and 0.8 (assuming ah ~ 800AU, 1,600AU and 

^ Tests found no significant difference betweeen using N=100 and 
N=300, so N=300 is assumed to be sufficient. 


3,OOOAU on the main-sequence), which is not a very strict 
constraint given the reduced perturbations due to the Galac¬ 
tic tide whilst the primary is on the main-sequence and the 
binary orbit tighter. Given the other uncertainities in our 
estimation of the fraction of white dwarfs polluted by the 
wide binary mechanism, we leave this constraint out of our 
estimation. 

4.2 The planetary system 

If /wB ~ 20% of wide binaries have orbital parameters such 
that a close approach occurs during the primary’s white 
dwarf evolution, the next question regards what fraction, 
fps, of these white dwarfs have planetary systems that can 
be scattered in such a manner as to produce white dwarf 
pollution. Our simulations (§3) necessarily consider a very 
specific architecture for the planetary system. Scattering, 
however, would occur in a wide range of planetary system ar¬ 
chitectures, including Oort cloud-like structures, rather than 
Kuiper-like belts. Clearly, the architecture of the planetary 
system has a great influence on the efficiency at which mate¬ 
rial is scattered inwards (e.g. Bonsor et al. 2012), however, 
for the purposes of this estimation we consider it sufficient 
to assume that all planetary systems where sufficient ma¬ 
terial survives have the potential to produce white dwarf 
pollution via the wide binary mechanism. An advantage 
of this mechanism is its lack of dependence on the exact 
structure of the planetary system. Outer planetary systems 
should survive the star’s evolution to the white dwarf phase, 
although planetesimal belts may become collisionally de¬ 
pleted (Bonsor & Wyatt 2010), and dynamical instabilities 
may clear some systems of material (e.g. Raymond et al. 
2009). Given the prevalence of planetary systems in the So¬ 
lar neighbourhood (e.g. Borucki et al. 2011; Mullally et al. 
2015) and that most planets or planetesimal belts cannot 
be detected, we consider a reasonable estimate to the frac¬ 
tion of wide binaries, where the white dwarf primary has 
a planetary system suitable for producing pollution to be 
fpg 50%. We note that this estimate is vague, and is ac¬ 
tually something that observations of polluted white dwarfs 
may help to better constrain. 

4.3 Conclusion 

We can estimate the fraction of white dwarfs with wide bi¬ 
nary companions likely to be polluted by this mechanism, 
at some point during their lives as white dwarfs, as 

fpoii ~ fwB X fps ~ 20% X 50% ~ 10% (5) 

with large uncertainties on our estimations of all of the pa¬ 
rameters. At any given instance in time pollution would be 
detectable in only a fraction, fr, of these systems. This frac¬ 
tion will depend on the timescale on which pollution is de¬ 
tectable following a close approach of the binary star. Many 
factors, of which are understanding is limited, need to be 
considered to produce a good estimate for this timescale, 
including the dynamical timescales in the planetary system, 
lifetime of any accretion disc formed, sinking timescale of 
metal pollutants, etc . We anticipate that the dynamical 
timescales in the planetary system will dominate, being sev¬ 
eral orders of magnitude longer than the other timescales 
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involved here. The dynamical timescales on which the outer 
planetary system continues to feed the white dwarf pollution 
may vary significantly depending on the exact architecture 
of individual planetary systems. Rather than attempt to ap¬ 
proximate these, we, therefore, conclude that a maximum of 
a few percent of any sample of white dwarfs with wide binary 
companions would be observed to have metal pollution pro¬ 
duced via the wide binary mechanism presented in this work. 
The wide binary fraction of white dwarfs is not well known, 
but if we consider that at least (and probably more than) 
10 % of white dwarfs should have wide binary companions, 
based on the 10% of solar-type stars observed with wide 
binaries (Raghavan et al. 2010; Lepine & Bongiorno 2007) 
(see discussion in the Introduction), then the wide binary 
pollution mechanism contributes to pollution in a maximum 
of a percent of any sample of white dwarfs. 


5 OBSERVATIONAL TESTS 

The mechanism presented here works only for a subset of 
the potential orbital parameters of the binary and plane¬ 
tary system. As discussed in §4, this fraction is likely to 
be of the order of a few percent of any random sample of 
white dwarfs with wide binary companions. Given that many 
observations find that more than 25%, and maybe up to 
50%, of white dwarfs are polluted (Zuckerman et al. 2003, 
2010; Koester et al. 2014), this wide binary mechanism is 
only ever going to make a minor contribution, and other 
mechanisms must be important in polluting white dwarfs. 
In fact, of roughly one hundred ‘pre-Sloan’ polluted white 
dwarfs, discovered by their proper motions, such that any 
companions is likely to have been found, only 5 have known 
common proper motion companions (Farihi et al, private 
communication). It is not clear whether this wide binary 
fraction is different from the general population, where, for 
example ~ 10% of solar-type stars are found to have wide 
separation (> 1,000AU) companions (Raghavan et al. 2010; 
Lepine & Bongiorno 2007). 

Zuckerman (2014) consider a sample of 17 white dwarfs 
with wide (> 2,500AU) companions, all observed with 
Keck/VLT to search for Ca II, and find that 5/17 (30%) are 
polluted. With such small number statistics, and significant 
bias in sample selection, it is difficult to make a good com¬ 
parison with observations of apparently single white dwarfs, 
where, for example 25% of DA white dwarfs exhibit Ca pol¬ 
lution (Zuckerman et al. 2003). It remains plausible that the 
wide binary mechanism could contribute to the pollution in 
these white dwarfs on the percent level, but without larger, 
well defined, samples, it is impossible to make a detailed 
comparison. 

The key evidence for a time-independent mechanism for 
pollution, such as suggested here, comes from the lack of 
any decrease in the fraction of polluted systems with age, 
and the observations of pollution in old (Gyr) white dwarfs, 
which, for example dynamical instabilities following stel¬ 
lar mass loss struggle to explain (e.g. Bonsor et al. 2011; 
Debes et al. 2012; Frewen & Hansen 2014). There are cur¬ 
rently insufficient white dwarfs with known binary compan¬ 
ions to assess any dependence in the fraction of polluted sys¬ 
tems with white dwarf cooling age. However, WD 1009-184, 
with its companion at 6,870AU and an effective tempera¬ 


ture of 9,940K (equivalent to a cooling age of ~700Myr) 
(Sion et al. 2009; Zuckerman 2014) stands out as a can¬ 
didate example of where the companion could potentially 
be responsible for the white dwarf pollution. Future obser¬ 
vational searches for companions to cool white dwarfs will 
provide critical evidence regarding the contribution of wide 
binaries to pollution. Gaia will play a crucial role in discov¬ 
ering companions to many nearby white dwarfs. 


6 DISCUSSION 

In this work we outline a proof of concept that illustrates 
the manner by which a wide binary companion could lead 
to pollution in white dwarfs. In §2 we illustrate that dur¬ 
ing the evolution of the binary’s orbit as it interacts with 
the Galactic tides, periods of increased eccentricity occur. 
These periods are most likely to occur during the white 
dwarf phase, due to the increase in separation of the bi¬ 
nary orbit following stellar mass loss. We claim that a wide 
binary companion whose eccentricity increases for the first 
time during the primary star’s white dwarf evolution has 
the potential to perturb any planetary system orbiting the 
white dwarf. In §3, we illustrated how planetesimals can be 
scattered onto star-grazing orbits, a pre-requisite for pollu¬ 
tion of the white dwarf’s atmosphere. This mechanism is 
independent of white dwarf cooling age, and could, there¬ 
fore, produce pollution around old (l-5Gyr) white dwarfs 
that other mechanisms struggle to explain. We estimate that 
this mechanism may contribute to the total fraction of pol¬ 
luted white dwarfs on the level of a few percent. There are 
no current observational samples that are sufficiently large 
to assess whether the population of polluted white dwarfs 
around single stars is different to those with companions. 

We consider that this mechanism is robust and will 
occur for white dwarfs with planetary systems and com¬ 
panions with the ‘correct’ parameters. The weakest part of 
this argument regards the link between planetesimals scat¬ 
tered inwards and the pollution of white dwarf atmospheres, 
which was not investigated in this work. The tidal disrup¬ 
tion of planetesimals and the acccretion of close-in dusty ma¬ 
terial onto white dwarfs have been previously investigated 
(e.g. Rafikov 2011a, b; Debes et al. 2012; Metzger et al. 
2012; Veras et al. 2014c, 2015c), although some gaps still 
remain in our understanding. Our discussion of whether a 
detectable level of pollution is produced is vague (see §4), 
and it remains possible that no detectable pollution is ever 
produced by planetesimal scattering. However, our simula¬ 
tions do show a significant increase in the number of plan¬ 
etesimals scattered onto star-grazing orbits with a binary 
companion on an eccentric orbit, compared to without a 
binary companion (see Fig. 5), and given the low accretion 
masses required to produce detectable pollution®, it remains 
plausible that the mechanism presented works for some re¬ 
gion of the parameter space. In fact, we anticipate that this 
mechanism has a weak dependence on the exact architecture 
of the planetary system, and a much stronger dependence on 
the orbital parameters of the binary. 

® Even one of the most highly polluted white dwarfs only needs 
a Ceres mass (1O“^M0) (Dufour et al. 2010), or see Girven et al. 
( 2012 ) 
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The purpose of this work was to outline a potential 
mechanism for the pollution of white dwarfs. We do not 
consider our estimation of the fraction of white dwarfs with 
wide binary companions polluted by this mechanism to be 
robust, and it is even trickier to determine the fraction of an 
observational sample that would be caught whilst displaying 
detectable pollution. Our estimation does, however, indicate 
that the wide binary mechanism is not the dominant mecha¬ 
nism for producing white dwarf pollution, even at late times, 
and that statistics in current observational samples are in¬ 
sufficient to determine whether or not this mechanism is pro¬ 
ducing pollution. Gaia and future observational searches for 
companions to nearby white dwarfs, should enable a clearer 
determination of its contribution. 

We have ignored the potential for exchange interactions, 
or changes to the binary’s orbit, which could have a big effect 
on the evolution of individual systems, but are not going to 
significantly affect the total population of white dwarfs with 
wide binary companions. We have also ignored any poten¬ 
tial evolution of the companion star. If it loses a significant 
fraction of its mass, this could increase the fraction of bi¬ 
nary orbits that meet our criteria for producing pollution, 
as in essence the system has double the chance of entering 
the correct regime. For example, one could imagine that the 
expansion in the orbit of the binary is insufficient to induce 
large perturbations from the Galactic tide as the primary 
loses mass to become a white dwarf, but becomes sufficiently 
wide once the secondary evolves to become a white dwarf. 
Perturbations could occur to planetary systems orbiting ei¬ 
ther the primary or the secondary. 

Our mechanism depends critically on the survival 
of a planetary system, containing sufficient material, to 
late times in the white dwarf phase. Although, theoreti¬ 
cally there are good arguments for the survival of outer 
planetary systems during post-main sequence evolution 
(Duncan & Lissauer 1998; Veras & Wyatt 2012), even if col- 
lisionally depleted (Bonsor & Wyatt 2010), there is a yet 
little direct observational evidence for the presence of plan¬ 
etary systems around white dwarfs (Mullally et al. 2008; 
Burleigh et al. 2008; Xu et al. 2015), although a few po¬ 
tential exceptions exist (Luhman et al. 2011; Marsh et al. 
2014). Planets have been detected around giant stars [e.g. 
Johnson et al. 2007, 2008; Sato et al. 2008) and horizontal 
branch stars (Silvotti et al. 2007). Dusty, debris-like discs 
have been observed at the centre of several planetary nebu¬ 
lae (Su et al. 2007; Chu et al. 2011), very early in the white 
dwarf phase. If future observations show that this mecha¬ 
nism is a significant contributor to white dwarf pollution, it 
provides further evidence for the survival of outer planetary 
systems to the white dwarf phase. 

The most important contribution of the wide binary 
mechanism presented here is to the pollution of old (Gyr) 
white dwarfs, that are hard to explain by other mecha¬ 
nisms. The mechanism is independent of white dwarf cool¬ 
ing age, except for any variation in the planetary system 
itself. On the main-sequence material can be depleted, both 
by collisions and dynamics, with a steep time dependence. 
By the white dwarf phase, collisional timescales are long 
(Bonsor & Wyatt 2010), and the system is much more likely 
to have reached a dynamically stable state, such that any 
such evolution is likely to only occur at low levels. 


7 CONCLUSIONS 

In this work we present a mechanism to explain the metal 
pollution observed in many white dwarfs, especially at late 
times. The orbits of wide binaries vary periodically with 
the Galactic tides. The influence of the Galactic tide in¬ 
creases following stellar mass loss. Thus, a wide binary may 
evolve for Gyrs at large separations, and yet, have a close 
approach for the first time during the primary’s white dwarf 
phase. Such a close approach could perturb any planetary 
system orbiting the primary, potentially scattering planetes- 
imals onto star-grazing orbits. 

We present simulations that illustrate the change in or¬ 
bital parameters of the binary (§2) and the scattering of 
planetesimals onto star-grazing orbits (§3). These demon¬ 
strations should be considered a proof of concept, rather 
than a detailed analysis. We estimate that this mechanism 
could contribute to pollution in up to a few percent of any 
sample of white dwarfs with wide binary companions, and 
up to a percent of all white dwarfs (§4), independent of the 
age of the white dwarf. Current observational samples do 
not have sufficient statistics to indicate whether or not this 
mechanism provides an important contribution to pollution 
in white dwarfs. The calcium-polluted, WD 1009-184, with 
its companion at 6,870AU and an effective temperature of 
9,940K (equivalent to a cooling age of ~700Myr) (Sion et al. 
2009; Zuckerman 2014) stands out as an example system 
where this mechanism could be acting. 

Wide binary pollution is unlikely to be the sole con¬ 
tributor to the population of metal-polluted white dwarfs. 
Instead, we envisage that it acts alongside other mecha¬ 
nisms, adding a time-independent pollution rate. Future ob¬ 
servational searches for wide binary companions to nearby 
white dwarfs {e.g. with Gaia) and an increased sample of 
metal-encriched old, white dwarfs are required to determine 
whether or not wide binary companions play an important 
role in white dwarf pollution. 
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